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Abstract 

In this paper we investigate relationships between the volumes of cells of three-dimensional 
Voronoi tessellations and the lengths and areas of sections obtained by intersecting the tessellation 
with a randomly oriented plane. Here, in order to obtain analytical results, Voronoi cells are 
approximated to spheres. First, the probability density function for the lengths of the radii of 
the sections is derived and it is shown that it is related to the Meijer G-function; its properties 
are discussed and comparisons are made with the numerical results. Next the probability density 
function for the areas of cross sections is computed and compared with the results of numerical 
simulations. 
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I. INTRODUCTION 



Three-dimensional Voronoi tessellations provide a powerful method of subdividing space 
in random partitions and have been used in a variety of fields such as computational geometry 



n °ri 



and numerical computing (see for instance [l| and the references therein), data analysis and 



n 



compression [2], geology |3|, and molecular biology [4, l5|. However, in many experimental 
conditions it is not possible to directly observe the cells themselves, just their planar or linear 
sections: thus it is of interest to study the relationships between the geometric properties of 



three-dimensional structures and their lower dimensional sections iGl. 

n 

A basic result has been proved in [7]: the intersection between an arbitrary but fixed 
plane and a spatial Voronoi tessellation is not necessarily a planar Voronoi tessellation. An 



analysis of some aspects of the sections of Voronoi diagrams can be found in jo]. However, a 
far as we know, no analytical formulas have been derived for the distributions of the lengths 
and areas of the planar sections: this is precisely the aim of this note. 

Following jsl, a three-dimensional Voronoi diagram will be denoted by V, and section of 
dimensionality s will be denoted by V(s,3). 

In the next section, some probability density functions used to fit empirical, or simulated, 
distributions of Voronoi cells size will be reviewed and in Section UTTj results will be presented 
on the distribution of the lengths and areas of V(2, 3). 

II. PROBABILITY DENSITY FUNCTIONS OF VORONOI TESSELLATIONS 

A Voronoi tessellation is said to be a Poisson Voronoi diagram (PVD), denoted by Vp, 
if the centers generating the cells are uniformly distributed. In the case of one-dimensional 
PVD, with average linear density A, a rigorous result can be derived, namely that the 
distribution of the lengths of the segments has the probability density function 

p{l) =AXHexp{-2Xl), (1) 
or, by using the standardized variable x = 1/ < I >= XI Q], 

= 4xexp (— 2x). (2) 
No analytical results are known for the distribution of the sizes of Voronoi diagrams in 
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2D or 3D; typically distributions of the surfaces or volumes of the Voronoi cells are fitted 
with a 3-parameter generalized gamma PDF of the standardized variable x, js], 10|. 



g{x; a, b, c) = c -""^ 



x"" exp (~bx'^ 



(3) 



r(a/c) 

It should be noted that from g{x; a, b, c) other, simpler, probability density functions can 
be derived, that have been also used to model empirical distributions of PVD: for instance, 
by setting a = b and c = 1, the one-parameter distribution proposed in js] results in 



h{x; b) 



xY ^exp^—bx), 



with variance 

Similarly, set a = b 
(131) then becomes 



m 



(4) 



(5) 



(3c? + l)/2, c = 1, where d is the dimensionality of the cells: Eq. 



p{x] d) 



3d+l 



a;~exp(-(3ci+ l)x/2). 



with variance 



0"„ 



3rf+ 1 



(6) 



(7) 



This PDF has been shown to give a good fit for Poisson Voronoi diagrams ll|], even 
though it has no free parameters: it will be used in the following to model distributions of 
the volumes of three-dimensional Voronoi tessellations. 



III. SECTIONS 

Consider a three-dimensional Poisson Voronoi diagram and suppose it intersects a ran- 
domly oriented plane 7: the resulting cross sections are polygons, as in the example shown 
in Fig. [H 

In order to obtain analytical results for the distributions of the lengths and areas of 
Vp(2, 3), some simplifying hypothesis is needed: here the cells will be supposed to be spheres. 
This is, admittedly, a rather rough approximation, but, on the other hand, it makes all cuts 
of the same simple shape, i.e. circles, irrespective of the orientation of the cutting plane. 
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FIG. 1: Planar cross sections of a three-dimensional Voronoi tessellation. 

As shown in the sequel, this approximation allows of deriving an analytical form for the 
distributions of the radii lengths and the areas of cell sections. 

Let fviy) be the probability density function for the cell volumes, then the PDF /ij(-R) 
for the lengths of their radii is given by 

fniR) = fv(^^rcR'^AnR\ (8) 

The use of Eq. ([H]) and of the standardized PDF with d = 3, leads to 

A plot of fn is shown in Fig. [2j 

Consider now the intersection of a sphere of radius R with a randomly oriented plane: 
the PDF of r, the length of the radius of the resulting circle, is jsj 

..(r) ^ f hm^^^dR (10) 
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FIG. 2: The PDF fa as a function of R. 

The probability of a plane's intersecting a sphere is proportional to i? [3| so, in conclusion, 
the PDF of r can be written as 



frir) = a / MR) 



-AR. 



(11) 



Insertion of Eq. (|9]) into (ITT!) results in the formula 

fr{r) = a 



4-10^ foo / 20 



243 



TT 



i?^*exp ( -—ttR^ 



-.dR. 



(12) 



The integral (fT2|) can be solved by making use of a computer algebra system such as MAPLE: 
the result is 

■5/6,1/6,1/2 \ ^^g^ 



Ur) = 2/3 a ^^^rGti(^ 



7/3,2/3,1/3,0,— ' ' 



where G is the Mejier G-function 12Nl4j: the numerical value of a, obtained via normal- 
ization of is a = 1.649. A plot of fr is shown in Fig. |3l The Meijer G function is very 
complicated and Eq. ( fT3l) does not lend itself to ready interpretation; moreover, at the best 
of our knowledge, it allows no simple approximations ( see 15|). 

However, some insight into the form of fr can be gained by considering Taylor expansions. 
In the interval [0, 0.1] fr is linear, fr{r) ~ 2.78r, while for r close to 1, in the interval [0.9, 1], 
it is well approximated by a quadratic polynomial: 



-0.006 (r-l) + 0.136 (r - 11 



(14) 
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FIG. 3: The probability density function fr as a function of r. 
The mode occurs in the interval [0.5, 0.6], where a close approximation of fr is given by 
fr{r) ^ 2.48 + 0.58 r - 95.528 (r - 0.55)^ - 202.18 (r - 0.55)''^ + 1901.87 (r - 0.55)^ ; (15) 

from this formula the mode m can be easily computed: m = 0.553. 

The mean and variance of fr are (r) = 0.487 and cr^ = 0.025, respectively; is moderately 
left skewed (skewness 7^ = —0.522) and is slightly platykurtic (excess kurtosis kr = —0.111). 
For comparison, note that fn is more symmetric ( skewness 7ij = 0.014) and has fatter tails 
(excess kurtosis kji = —0.005). 

The distribution function Fr is 

F(r)-—a 35/61 02/3^4,2 f 100 o fi 1,7/6,1/2,5/6 \ J_ gx 
~ 90 V 9 8/3,i,2/3,i/3,f ,0 j y^^) 

The PDF Ja of the areas of Vp(2, 3) can be obtained from fr by means of the transfor- 
mation 



= M ( 4V 1 ^^-'^ (17) 



TT 



that is, 



/.(^)= 0,549^^*^G«(H^|V../«/.^„),-^ . (18) 

Since, for r close to 0, fr{r) ~ r from Eq. ( fTTl) it follows that /a(0) 7^ and this result 
has been also verified by numerical simulations we have carried out (see also [6] for a further 
example). Indeed it is easy to verify that /yi(0) = 0.443. Fig. H] shows the graph of Ja- 
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FIG. 4: The PDF Ja as a function of A. 

The mode of Ja is mA = 0.858, its mean and variance are (A) = 0.824 and = 0.204, 
respectively; is approximately symmetric (skewness 7^ = 0.278) and its excess kurtosis 
is kA = —0.337, indicating that the transformation from fr to Ja yields a PDF with a lower 
and broader peak and with shorter and thinner tails. 

The distribution function Fa is given by 

9 vr ' ' ! ~ 

A comparison has been carried out between the distribution functions Fr, Fa, as as given 
by Eq. (|T6l) . ( fT9l) respectively, and the results of a numerical simulation, performed as 
follows: starting from 300000 3-dimensional cells 100168 irregular polygons were obtained 
by adding together results of cuts by 41 triples of mutually perpendicular planes. 

The area of the irregular polygons has been obtained with the algorithm described in 



FA = 0.0l83'/H0'/'G'Ai— — 



161] : briefly, each cut defines a two-dimensional grid and the polygons areas we measured 
by counting the number of points of the grid belonging to a given cell. 

As concerns the linear dimension, in our approximation the two-dimensional cells were 
considered circles and thus, for consistency, the radius r of an irregular polygon was defined 
as 

r^(-) . (20) 

that is r is the radius of a circle with the same area of the polygon, A. The empirical 




FIG. 5: Comparison between theoretical distribution (continuous line) and data (empty circles) 
for the distribution of radii r of 2D cells which result from an intersection with a plane. The 
maximum distance between the two curves is dmax = 0.044 . 

distribution of radii is shown in Fig. |5] together with the graph of F^: the maximum distance 
between the empirical and analytical curves is dmax = 0.044. Likewise a comparison between 
Fa and the simulation is shown in Fig. |6l in this case dmax = 0.039. 

IV. CONCLUSIONS 

In this paper analytical formulas have been provided which model the distributions of 
the lengths and areas of the planar sections of three-dimensional Poisson Voronoi diagrams: 
in particular, it has been shown that they are related to the Meier G function. 

This finding is consistent with the analytical results presented in |l17|], where it is proved 
that nonlinear combinations of gamma variables, such as products or quotients, have distri- 
butions proportional, or closely related, to the Meijer G distribution. 

The analytical distributions F,. and F^ been compared with results of numerical simula- 
tions: in evaluating the differences between analytical and empirical distributions it must 
be kept in mind that the cells were approximated as spheres and that the distributions 
used here have no free parameters that can be adjusted to optimize the fit. The results 
obtained here may be useful for applications in stereology in that they allow of predicting 
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FIG. 6: Comparison between data (empty circles) and theoretical curve (continuous line) of the 
distribution of areas. The maximum distance between the two curves is dmax = 0.039. 

the distributions of linear and planar measures of sections, given an arrangement of three- 
dimensional cells. This method can also have application to astrophysics, namely in the 
analysis of the spatial distributions of the voids between galaxies, that are experimentally 



measured as two-dimensional slices 18| and that can be identified with two-dimensional 
cuts of a three-dimensional structure. For instance, the distribution of the radius between 
galaxies of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) has been reported 
in jl8|; this catalog contains 1054 radii Vg of voids. To compare these experimental data 
with the results obtained here one should know the values of Re, the radii of voids in three 
dimensions; however by making use of standardized variables x = r/{r), Xe = Te/(?"e) one 
obtains for the standard deviation based on observations a^^ = 0.26, close to the theoretical 
value, = 0.32. Finally, the 2D cuts presented in this paper may be related to Laguerre de- 



compositions (weighted Voronoi decomposition) |l9 



42 1|. with a weighting factor depending 



on the distance between the between cells and the cutting plane. 
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